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ASSIGNING FREQUENCY CHANNELS IN A CELLULAR SYSTEM 

This application claims the benefit of U.S. Provisional Application No. 

60/ , , titled "Assigning Frequency Channels In A Cellular System", by Krishnan 

Kumaran and Boris D. Lubachevsky, filed May 10, 2001. 

BACKGROUND OF THE INVENTION 

Field of the Invention 

This invention relates generally to wireless communications and, more 
particularly to cellular wireless systems. 

Discussion of the Related Art 

In a cellular system, a coverage region includes a plurality of coverage cells in 
which different base stations service cellular communications. Each base station has a set 
of frequency channels to support communications with mobile units located in the one or 
more cells serviced by the station. Since a cellular system licenses a fixed frequency 
spectrum, the channels of the spectrum must be shared by different users. 

Different ways of sharing frequency channels among the various base stations of 
the system often provide different call and data capacities. To increase the call and data 
capacity, the same frequency channel is available to support more than one simultaneous 
cellular communication. The level of such channel reuse is however, limited by 
interference that the reuse causes between different communications. Interference occurs 
among different communications supported by the same frequency channel, e.g., to 
adjacent base stations, and among different communications supported on adjacent 
frequency channels. One goal in the design of cellular systems is to increase system 
capacities by enabling more reuse of frequency channels by base stations while 
minimizing interference between independent communications. 

Some cellular systems implement time division multiple access (TDMA) to 
support communications between base stations and cellular users. In TDMA, a base 
station defines a series of time slots in individual available frequency channels. The base 
station uses the different time slots in a frequency channel to support simultaneous 
communications with different users. 
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BRIEF SUMMARY OF THE INVENTION 

Various embodiments prioritize the frequency channels of a cellular system for 
servicing communications based on learning. In particular embodiments, the learning 
results from processing communication traffic by a dynamical algorithm that evolves lists 
of channel rankings associated with base stations in a manner that reduces overall inter- 
communication interference. The lists prioritize the channels for use in servicing 
communications by the base stations. The algorithm evolves the channel rankings 
towards a stable fixed point. The fixed point provides lists of channel rankings that 
enable the same channels to be used by different ones of the angular sectors without 
above-threshold interference between independent communications. The fixed-point 
values of the channel rankings provide for better operation of the cellular system, i.e., 
high capacity and/or lower interference, than the initial channel ranges and typically, are 
not strongly dependent on the initial channel rankings. 

In one aspect, the invention features a process for assigning channels to service 
communications in a cellular wireless system. The process includes performing a 
simulation of the system to produce a plurality of lists of channel rankings and assigning 
the lists of produced channel rankings to the base stations for use in assigning channels to 
service communications with mobile units. The simulation evolves the lists of channel 
rankings according to an algorithm that dynamically reduces inter-communication 
interference. 

In another aspect, the invention features a program storage medium encoded with 
a computer executable program that is configured to perform the above-described 
process. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 shows an exemplary cellular wireless system; 

Figure 2 shows a list that ranks frequency channels for servicing communications 
by one angular sector of Figure 1 ; 

Figure 3 schematically illustrates a system that determines rankings of frequency 
channels in the cellular wireless system of Figure 1; 

Figure 4 is a flow chart for a process that provides lists that rank frequency 
channels for use by the angular sectors in the cellular system of Figure 1 ; 
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Figure 5 is a flow chart illustrating a program that simulates the evolution of lists 
of frequency channel rankings dynamically; 

Figure 6 shows data structures used by the program of Figure 5 to dynamically 
evolve lists of frequency channel rankings for the cellular wireless system of Figure 1; 

Figure 7 is a flow chart that shows how the program of Figure 5 processes a new 
communication creation event; 

Figure 8 is a flow chart that shows how the program ov Figure 5 processes a 
communication dialing or set up event; 

Figure 9 is a flow chart that shows how the program of Figure 5 processes a 
communication maintenance event; 

Figure 10 is a flow chart that shows how the program of Figure 5 processes a 
communication termination event; 

Figure 1 1 is a flow chart that shows how the program of Figure 5 processes an 
event to re-rank or reprioritize usage of frequency channels; 

Figure 12 is a flow chart that shows how the program of Figure 5 processes an 
event to determine the stability for the cellular system's channel ranking lists; and 

Figure 13 is a graphical representation of the evolution of a quantity whose values 
are indicative of convergence to a fixed point during simulation of a cellular system. 

DETAILED DESCRIPTION OF THE EMBODIMENTS 

A dynamical system can incorporate learning into a system's temporal evolution. 
The various embodiments dynamically evolve a cellular system according to a dynamics 
to rank channels for use in servicing communications, i.e., calls to and from mobile units, 
in a manner that produces less interference. The dynamics usually has one or more fixed 
points that provide lists of channel rankings for use by the actual cellular system. 

Figure 1 shows a portion of a cellular wireless system 2. The cellular wireless 
system 2 includes mobile units Mi-M m and wireless base stations Bi-B„ that are able to 
communicate with nearby ones of the mobile units Mi-M m . The base stations Bi-B n have 
fixed locations that are spread over the system's coverage region. Each base station Br 
B n has a set of fixed angular sectors Sj-Sk in which the base station is able to maintain 
independent communications with mobile units Mi-M m . The mobile units Mi-M m are 
located and may move in the coverage region. At a particular time, one of the mobile 
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units Mi-M m may or may not have active communication set up with one of the base 
stations Bi-B n . 

Each angular sector Si-Sk of a base station Bi-B„ describes a range of directions in 
which the particular base station Bi-B n is able to maintain independent communications. 
For example, base station Bi has angular sectors defined by: 0 = 0 0 - 120 °, 0 = 120 0 - 
240 °, and 0 = 240 0 - 360 °. In exemplary embodiments, the number, angular size, and 
orientation of the various angular sectors may differ for different base stations Bi-B n . 

Each angular sector Si-Sk has an associated transmission power, reception 
sensitivity, and a list of channel usage rankings. The power available for 
communications with a particular mobile unit Mi-M m depends on the sector's 
transmission power and a fading matrix. Similarly, the sensitivity available for 
communications with a particular mobile unit Mi-M m depends on the sector's reception 
sensitivity and the fading matrix. The fading matrix accounts for signal attenuation 
properties between a particular base station Bi-B n and a particular location, i.e., the 
location of a mobile unit Mi-M m . The power, sensitivity, and channel ranking may differ 
for the angular sectors Si-Sk of the same base station Bi-B n . 

Different ones of the angular sectors Si-Sk reuse the same set of frequency 
channels. While reuse of the same or adjacent frequency channel leads to interference 
between different cellular calls, channel reuse also increases the traffic bandwidth of the 
cellular system 2. The amount of inter-communication interference depends on 
proximity's of and transmitted power in angular sectors Si-Sk that use the same and/or 
adjacent frequency channels during the same time slots. 

For one angular sector Si-Sk, a channel ranking is a list that prioritizes the 
frequency channels of cellular system 2 by their quality for supporting communications 
between the sector and a mobile unit Mi-M m . Figure 2 shows a list 3 for an exemplary 
channel ranking of one angular sector. The list 3 has entries 4-7 that provide identifiers 8 
and quantities, |i n , representative of measured interference levels 9 of the frequency 
channels for communications serviced by the particular angular sector. Typically, the 
j^/s are only "representative" of interference levels, because they are defined in a manner 
that smoothes out rapid fluctuations in interference caused by the usage of frequency 
channels (see below). The entries 4-7 are ordered so that the channel with the lowest 
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inter-communication interference level, i.e., lowest is first and the channel with the 
highest inter-communication interference, i.e., highest \in, is last. 

Channel usage by an associated angular sector is controlled by the channel 
rankings of list 3. The rankings enable the associated base station Bj-B n to find an 
available channel with lowest interference by searching from the top of the list 3, i.e., the 
rankings prioritize the channels for use in supporting communications. Such a search can 
also be used to determine whether the best available channel provides too much 
interference. Thus, the search determines whether the channel has an above-threshold 
interference. 

For cellular system 2, a complete channel ranking includes "k" instances of lists 
of channel rankings. Each instance is a list similar in form to the list 3 of Figure 2. Each 
list ranks the system's frequency channels for one of angular sectors SpSk of therein. In 
various embodiments, the complete set of lists of channel rankings is dynamically 
evolved to find a set of lists of rankings that produce low interference for the whole 
cellular wireless system 2. 

Figure 3 schematically illustrates a system 12 that determines and updates the lists 
of channel rankings for cellular wireless system 2 of Figure 1. The system 12 receives 
input data from base stations Bi-B n and sends new channel assignments to the base 
stations Bi-B n via link 18. The input data characterizes relevant aspects of the system 12 
for determination of the rankings. The system 12 supplies new lists of channel rankings 
that compensate for changes to system properties such as communication traffic levels 
and mobile unit distributions. A base station Bi-B n assigns a new communication to the 
one of its highest ranked channel(s) angular sector Si-Sk that is best oriented to support 
the communication. 

The system 12 includes a microprocessor 14 with an internal data storage device 
15, e.g., a hard drive or an active memory, and an external drive 16. The microprocessor 
14 executes a computer program that is stored either in internal device 15 or in an 
external data storage media 17, e.g., an optical or magnetic disk that is executable in 
drive 16. The program determines a list of channel rankings for each angular sectors Sj- 
Sk of cellular wireless system 2 of Figure 1. After determining the lists of channel 
rankings, the processor 14 sends the lists to the base stations BpB„ via link 18. Each 
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base station Bi-B n stores the lists for the channel rankings for its own angular sectors Si- 
S k on an internal data storage device MBi-MB n . To change the channel rankings, e.g., in 
response to the receipt of new input data from the base stations Bi-B n , the system 12 
sends new lists of channel rankings to one or more of the base stations Bi-B n . 
5 Figure 4 is a flow chart for process 20 performed by a channel allocation system, 

e.g., system 12 of Figure 3, to find lists of channel rankings for use in cellular system 2 of 
Figure 1 . The allocation system receives input data characterizing the cellular system 2 
for which the channel rankings will be made, e.g., data characterizing a city and its 
suburbs (step 22). The input data includes the cellular system's frequency channel set; 
10 the base stations locations; the orientations of angular sectors of the base stations; the 
local densities and communication traffic rates for mobile units; a fading matrix; and the 
communication quality requirements. The fading matrix describes signal attenuation 
values between pairs of base stations and mobile unit locations. The communication 

^ quality requirements fix upper bounds on acceptable levels of inter-communication 

few 

*S 15 interference. 
Q 

p An exemplary channel ranking process uses uplink communications, downlink 

J^j communications, or both to determine the lists of channel rankings. For uplink and 

*p downlink based systems, the input data also includes reception sensitivities of angular 

sectors and transmission powers of the various angular sectors, respectively. 
20 The input data may also include reception sensitivities and/or transmission powers 

^ of the mobile units. 

J^f In some embodiments, additional input data is received by the system that 

determines the channel rankings. One such type of data characterizes time division 
multiplexing (TDM) performed by the angular sectors of cellular system 2. For example, 

2 5 TDM data indicates relative time shifts between the start times of multiplexing time slots 

of the different angular sectors. Another such type of data includes a set of initial lists of 
channel assignments. For example, this type of data is provided to system 12 of Figure 3, 
because the system 12 updates the lists of channel rankings to follow time-dependent 
changes in the cellular system 2. Another such type of data characterizes durations of 

3 0 communications. Another such type of data characterizes redial rates for block 

communications. Another such type of data characterizes interference levels or power 
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levels that cause handoff initiations. Another such type of data may describe motion of 
mobile units Mi-M m through the time evolution of local densities for the mobile units 
Mi-M m . 

In some embodiments, the input data provides for several classes of mobile units 
5 that are distinguished by communication quality requirements, motion, redial rates, 
communication durations, interference tolerances, and/or communication bandwidth 
requirements. 

From the input data, the channel allocation system siriiulates cellular system 2 of 
Figure 1 by a dynamics that evolves the frequency channel rankings to reduce inter- 
10 communication interference, i.e., a dynamics that learns how to reduce interference (step 
24). The dynamics evolves "k" lists that rank the same set of system frequency channels, 
e.g., lists like list 3 of Figure 2, with respect to their usage by the various angular sectors 
Si-Sk. At any particular simulation time, the channel allocation system stores the lists of 
channel rankings, which have been produced by the simulated dynamics for each angular 

*S 15 sector S i -Sk of the cellular system (step 26). 

O 

The simulation evolves channel rankings to reduce overall inter-communication 
interference while maintaining call blocking and call drop rates below selected threshold 
values. Thus, the stored lists of channel rankings typically, provide lower interference 
L» levels than the initial lists without producing high blocking or drop rates. The channel 

M 2 0. allocation system sends the stored lists of channel rankings for the angular sectors Si-Sr 

□ 

jrj to the base station Bi-B n on which the sectors Si-Sk are located, e.g., via link 1 8 in Figure 

2 3 (step 28). 

The received channel rankings, which are lists similar to list 3 of Figure 2, control 
subsequent channel usage decisions by base stations Bi-B n . In particular, the base 

2 5 stations Bi-B n assign new communications to or from mobile units Mi-M m to unused 

communication time-slots of the highest ranked channels based on the rankings of the 
lists associated with the angular sectors Si-Sk used to support the communications. 

Figure 5 is a flow chart illustrating one program 30 that implements the simulated 
dynamics of Figure 4. First, the program receives a set of initial lists of channel rankings 

3 0 for the angular sectors Si-Sk (step 32). The initial rankings include "k" lists of channel 

rankings. Each list ranks the same system frequency channels according quantities, |in, 
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that are representative of interference that an associated angular sector Si-Sk would 
encounter in using the channels to support communications. The initial rankings are 
fixed either randomly or based on input data. For random initial rankings, interference 
levels are initially assigned in a manner that is consistent with the channel assignments, 
i.e., higher ranked channels are assigned lower values for the quantities, ^i n , that 
characterizes inter-communication interference. 

After the initial assignments, the program 30 evolves the channel assignments 
dynamically in a manner that lowers inter-communication interference (step 34). The 
program processes simulated mobile communications to establish populations of used 
frequency channels, i.e., channels supporting calls with mobile units, for each angular 
sector Si -Sic of the system 2 being simulated. The used channels and used time slots 
thereof for the various angular sectors Si-Sk characterize a time-dependent system state. 
From the system-state, the program 30 measures interference levels, i.e., values of m n , 
that would occur between communications over various channels for the different angular 
sectors Si-Sk. Based on the measured interference levels, the program 30 modifies 
channel rankings, which modifies, in turn, subsequent channel usage by the angular 
sectors Si-Sk in the simulation. Thus, the program 30 "learns" to reduce interference 
dynamically while processing simulated communications with mobile units Mi-M m . 

The program also regularly makes measurements to determine whether the 
simulated evolution has neared a stable fixed point from which further evolution of the 
lists of channel rankings does not occur. In response to being near a fixed point, the 
program stores the instantaneous lists of channel rankings for the angular sectors Si-Sk to 
a separate data storage file from which they can be read out (step 26'). The channel 
allocation system sends these lists of channel rankings to the base stations Bt-B n at step 
28 of Figure 4. 

Figure 6 shows data structures used by one embodiment of program 30 uses to 
perform an event-driven dynamics. The data structures include an event queue 38 and a 
system-state file 39. The event queue 38 is a sequence of events that controls the 
evolution of a discrete-time dynamics simulated by the program. The event queue 38 
includes events related to communications, e.g., new call creation, call dialing, call 
redialing, call maintenance, call termination, call handofF, and events related to 
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dynamical evolution of the lists of frequency channel assignments, e.g., channel re- 
ranking and stability determination. The system-state file 39 stores data indicating 
presently used channels and time slots for angular sectors Si-Sk, densities of available 
mobile units Mi-M m , and other data defining the instantaneous state of the simulated 
system as described below. 

The program 30 provides a simulation dynamics that is controlled by event queue 
38. The dynamics involves sequentially processing events in the queue 38 and updating 
the system-state file 39 and queue 3o during the event processing. At any time, a 
program pointer 40 indicates the event presently being processed, e.g., event K in Figure 
6. In the dynamics, new call creation, channel re-ranking, and stability determination 
events are either initially added to the event queue 38 or periodically added to the event 
queue 38. Processing of these three types of events results in the insertion of new events 
into the event queue 38. Termination of the simulation occurs when the processing of a 
stability determination event indicates that the evolution of the lists of channel rankings 
has approached close to a stable fixed point of the evolution. In response to such a 
determination, the program stores the lists of channel rankings as described in step 26' of 
Figure 5. 

The dynamics that program 30 uses to process events related to calls and events 
related to channel rankings are described by flow charts shown in Figures 7-10 and 
Figure 11-12, respectively. The process of the various event types includes inserting, 
modifying, and/or removing later events from event queue 38 and/or updating the 
system-state file 39, e.g., changing designations of the used frequency channels and times 
slots therein for the system's angular sectors. 

Figure 7 is a flow chart 50 showing how program 30 processes a new call creation 
event. In response to the new call creation event, the program 30 selects a delay time for 
an attempt to set up the new call, i.e., to dial the new call (step 52). The inter-arrival 
times "T" between the future times of the dialing event are distributed according to a 
probability distribution, P(T), e.g., a Poisson distribution. The distribution may be part of 
the input data or be a default distribution, e.g., P(T) = [a]exp(-aT). In response to the new 
call creation event, the program 30 also selects a location for the mobile unit Mi-M m that 
will participate in the new simulated call (step 54). The mobile unit's location is 
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determined from local densities, and call traffic rates of the mobile units Mi-M m . In 
some embodiments, the determinations of mobile unit locations also depends on 
availabilities, i.e., mobile units involved in other calls are unavailable. In the simulated 
dynamics, these properties depend both on the original input data and on the system's 
state as defined by system-state file 39. For example, the availability of a mobile unit 
M r M m depends on whether the unit has a call, and the local density of mobile units M r 
M m depends on mobile unit motion in embodiments in which the mobile units Mi-M m 
move. 

Some embodiments use stationary densities of mobile units Mi-M m to perform the 
simulation based on call creation events that fix the call dialing time and mobile unit's 
location. A new event to dial the call is inserted into event queue 38 at an entry location 
that corresponds to the call dialing time (step 56). The new event indicates both the 
location of the mobile unit Mi-M m that will dial or receive the call and the call dialing 
time, i.e., position in the queue 38. 

Figure 8 is a flow chart 60 showing how program 30 of Figure 5 processes a call 
dialing event, i.e., an event to set up a scheduled call. In response to such an event, the 
program 30 ranks angular sectors Si-Sk of base stations Bi-B n located within a 
predetermined distance of the mobile location (step 62). The angular sectors Si-Sk are 
ranked based either on a sectors 1 abilities to deliver signal power to the location of the 
mobile unit Mi-M m participating in the call or on sectors 1 sensitivities to signal power 
from the location of the mobile unit Mi-M m . For an angular sector Si-Sk, the sector's 
ability to deliver signal power to or sensitivity to signal power from a location is 
determined from input data that fixes the sector's power or sensitivity and from a fading 
matrix. 

The fading matrix describes signal attenuation between particular base stations 
Bi-B n (row index) and particular locations (column index) of the mobile units Mi-M m . 
Often, the fading matrix is based on a power law attenuation of transmitted signal power. 
An exemplary law provides that transmitted signal power is attenuated like (distance)" 3 5 . 
The entries of the fading matrix also account for shadow attenuation, e.g., as caused by 
buildings or mountains. In exemplary embodiments, the entries of the fading matrix are 
derived by program 30 from the input data. 
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Program 30 selects the highest ranked angular sector Si-S k that is available to set 
up the call being dialed (step 64). Then, the program 30 determines whether the selected 
angular sector Si-Sk has one or more acceptable time slots in a frequency channel to 
support the dialed call with the participating mobile unit (step 66). For the selected 
angular sector Si-Sk, the program 30 searches for unused communication time slots in the 
frequency channel with the highest ranking, i.e., the channel with the lowest interference 
level. The highest ranked channel is found from the system-state as defined by system 
state file 39. 

A potential time slot is only acceptable if the ratio of the communication power in 
the slot to the interference power in the same slot (i.e., in the same and adjacent 
frequency channels) is greater than a predetermined threshold. The communication 
power depends on the power of the participating mobile unit or angular sector and on the 
appropriate element in the fading matrix. 

In searching for unused time slots, the program 30 considers up to a preselected 
number of the highest ranked channels in the selected angular sector Sj-Sk. If a time slot 
is found for the call, the program 30 assigns the time slot to the new call by updating the 
status of the time slot in the system-state file 39 (step 68). In this case, the program 30 
finishes setting up the call by creating a future call maintenance event and a future call 
termination event in event queue 38 (step 70). When processing a call maintenance 
event, a decision event is made as to whether the call in progress should be continued, 
dropped, or handed off. Drops and hand-offs are based on call quality considerations. 
For example, a call drop occurs if the ratio of the signal power to the interference noise 
power too large at the time of the maintenance event. Similarly, in embodiments 
implementing handoffs, a handoff occurs when the call quality is improved by 
transferring the call to another base station Bi-B n . The call termination event causes 
termination of the call after a selected time, e.g., defined by a probability distribution 
whose mean time is part of the input data. 

If a time slot is not available in the selected angular sector Sj-Sk, program 30 
determines whether other angular sectors Si-Sk remain in the set of ranked angular 
sectors (step 72). If other angular sectors remain, the program 30 loops back 74 to repeat 
the process for the highest ranked one of the remaining angular sectors Si-Sk. If no 
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ranked angular sectors Si-Sk remain, the program 30 creates in event queue 38 an event to 
redial the call at a future time, because the call is temporarily blocked due to the absence 
of an appropriate angular sector to support the call (step 76). The redial event simulates a 
caller's persistence to set up the new call, i.e., by an attempt to redial the call at a future 
time. 

Program 30 processes redial events like new call dialing events except that redial 
events have associated persistency ratings. The persistency rating tracks the number of 
tries to redial the blocked call. Afc^r a selected number of redial tries, the program 30 
does not create further call redialing events. The persistency rating simulates a caller's 
persistence in retrying to connect a blocked call. In some embodiments, redial 
persistencies are part of the input data. 

Figure 9 is a flow chart 80 showing how program 30 processes a call maintenance 
event. In response to such an event 90, the program 30 determines whether the quality of 
the communication channel that is supporting the call is above a predetermined threshold 
value (step 82). The quality is determined from a ratio of the signal power delivered by 
the servicing angular sector Si-Sk to the participating mobile unit Mi-M m , or the signal 
power received from the participating mobile unit Mi-M m at the servicing angular sector 
Si-Sk to the interference power in the channel of the angular sector Si-Sk supporting the 
call. The power levels are based on input data on properties of the angular sector SpSk 
that is servicing the call, properties of the mobile units Mi-M m , and the fading matrix. 
The power levels can change, e.g., due to simulated motion of the mobile units Mi-M m . 
If the ratio of the relevant signal power to the interference (i.e., noise) power level is 
above a predetermined threshold, the program 30 maintains the call and creates a future 
call maintenance event for the call (step 84). If the ratio of the relevant signal power 
level to the interference power level is below the predetermined threshold, the program 
30 marks the call as bad and tries to hand off the call to a different angular sector/base 
station if the simulation supports handoffs (step 86). If the hand off fails after a selected 
number of the tries the call is dropped. If the embodiment does not support handoff, the 
call is simply dropped. 

Embodiments that support handoffs treat a handoff situation like an event to dial a 
new call, i.e., to set up the call with a new angular sector Sj-Sk. Handoffs are, e.g., 
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supported in embodiments of program 30 that simulate motion of mobile units Mi-M m . If 
a handoff is not possible, the call is terminated. A termination removes other events 
related to the same call from event queue 38, e.g., based on a listing of the other events 
that is stored in system-state file 39, and updates the system-state file 39 to show that 
time slots that previously supported the call are now unused. 

Figure 10 is a flow chart 90 showing how program 30 processes a call termination 
event. In response to such an event, the program 30 updates syst-sm-state file 39 to free 
the time slot previously used to support the call (step 92). Freeing the time slot involves 
updating system-state file 39 to show that the time-slot of the frequency channel in the 
angular sector Si-Sk that previously supported the call is now unused. The program 30 
also removes from event queue 38 any future events that are related to the same call, i.e., 
based on a list of such events that is stored either in the system-state file 39 or in the 
event queue 38 itself (step 94). In some embodiments, the program 30 also updates the 
system-state file 39 to show that the mobile unit Mi-M m previously participating in the 
call is now available for other calls, i.e., the program 30 updates the density of available 
mobile units Mi-M m . 

Figure 1 1 is a flow chart showing how program 30 processes an event to re-rank 
one or more frequency channels of an angular sector Si-Sk* The program 30 selects an 
angular sector Si-Sk whose ranking is to be updated (step 102). In exemplary 
embodiments, this selection of angular sectors Si-Sk is either performed randomly or in 
round robin fashion. For the selected angular sector Si -Sk, the program 30 recalculates a 
quantity, \x n , representative of interference that a call supported by one or more selected 
frequency channels of that selected angular sector Si-Sk would experience (step 104). 
The program 30 updates the system-state file 39 with the recalculated |On values (step 
106). 

The calculation of |Xn at time "n" is based on a measured interference level, i.e., 
m n in decibels, and on the previous value of jOn-i for the same frequency channel of the 
angular sector under consideration. The measured interference level, m n , depends on 
which other channels in the selected angular sector and on which channels in other 
angular sectors are presently supporting communications. Typically, interference 
measurements are restricted to nearby base stations Bi-B n and to simulated 
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communications in the same and adjacent frequency ranges, because other 
communications produce little interference for the selected angular sector and channel 
thereof. Interference levels are averages over various time slots of the selected channel in 
simulations that simulate time division multiplexing. The contribution to ja n from the 
previous interference value, i.e., |a n -h smoothes out rapid fluctuations in channel activity 
levels. In exemplary embodiments, ^ n = (1- ot)|i n -i + cmi n with 0 < a < 1 and preferably 
<x<0.1. 

From the updated list of values of |i m the program 30 updates the rankings of 
frequency channels for the selected angular sector Si-Sk (step 108). This updating 
includes updating the list of channel rankings that is associated with the selected angular 
sector Si-Sk. The lists of channel rankings for the various angular sectors Si-Sk are 
located in the system-state file 39, e.g., list 3 of Figure 2. The updated list of channel 
rankings affects the future processing of call dialing events. 

In various embodiments, program 30 re-ranks different portions of the list of 
frequency channels in response to an event to re-rank the channels for a particular angular 
sector. In some embodiments, the list of channel rankings for each angular sector 
includes a pointer so that updating of the \i n values is performed in round-robin fashion 
for individual sectors. Then, processing of several events to re-rank the channels of the 
same angular sector results in an update of the entire list of channel rankings of the 
sector. 

In some embodiments, each list of frequency channels, e.g., list 3 of Figure 2, 
includes a short list of the highest ranked channels for the associated angular sector Si-Sk. 
The short list includes the subset of the channels that the associated angular sector Si-Sk 
uses to support communications with mobile units Mi-M m . Typically, events to re-rank 
frequency channels more frequently re-rank the channels on the short list, channels and 
less frequently re-rank the remaining channels. In one embodiment, the (in's that are 
indicative of interference levels evolve differently for channels on the short list and for 
the remaining channels. In particular, the equations, i.e., |Xn = (1- a)|0n.i + ct|in> use an a 
= as in the equations for channels on the short list and use an a = <xl in the equations for 
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the remaining channels. Furthermore, ct s > > a L , which implies that the channels on the 
short list have their |i n variables updated more rapidly. 

Figure 12 is a flow chart 1 10 showing how program 30 processes a stability 
determination event, which is indexed by a temporal sequence number "q". In response 
to such an event, the program 30 evaluates a quantity, Rq whose value is indicative of 
whether the simulation of the system f s channel rankings is converging to a fixed point 
(step 112). An exemplary formula for Rq is: 

Rq = X angular sectors Sp (rank change of last updated channel)/k. 
The rank changes are integers that describe how far the last updated frequency channel's 
rank was changed during the last re-ranking event for the associated angular sector Si-Sk. 
The processing of a channel re-ranking event includes storing this integer in a subfile 
located in the system-state file 39. The subfile has k entries, i.e., one entry for each 
angular sector Si -Sr. 

Figure 13 shows the general behavior of the quantity Rq with the event index 
number "q". Initially, Rq is large and may increase or strongly fluctuate with increasing 
sequence number "q". But, as a fixed point is approached, fluctuations in the value of Rq 
start to decrease, because the re-ranking process of Figure 1 1 leads to less reorderings of 
the frequency channels for the various angular sectors SpSk. As the simulation time 
increases, the values of Rq for sequential "stability determination events", i.e., sequential 
"q", go to a constant value, e.g., zero. The values of the Rq's also become smaller as "q" 
grows, because fewer re-ranking of frequency channels occur. Of course, other formula 
for quantity Rq could also provide measures of the distance of the system evolution from 
a fixed point. 

Referring again to Figure 12, program 30 compares the value of quantity Rq for 
the latest stability determination event "q" with the value of Rq. p for an earlier stability 
determination event "q-p" where "p" is a preselected number (step 1 14). The comparison 
determines whether |Rq - Rq. p | is converging to zero. In some embodiments, the program 
30 finds convergence if |Rq - Rq. p | has monotonic decreasing oscillations as a function of 
"<T> e.g., oscillations whose magnitude decrease faster than a preselected function of "q". 
In other embodiments, the program 30 finds convergence if |Rq - Rq. p | is smaller than a 
predetermined threshold. In response to finding such convergence, the program 30 stops 
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the simulation and records a list of the channel rankings for each angular sector Si-Sk 
(step 116). In such a case, the dynamical simulation is near a fixed point for dynamical 
evolution of the "k" lists of channel rankings. Then, the recorded lists of channel 
rankings provide data that base stations Bi-B n of the actual cellular system 2 will use to 
assign frequency channels during later processing of real communications with mobile 
units Mi-M m . If the program 30 does not find that |Rq - Rq. p | is converging, the program 
30 continues the simulation of cellular system 2 (step 118). 

Other embodiments of simulation program 30 handle events related to call 
processing and events related to the evolution of the lists of channel rankings separately. 
For example, one embodiment includes separate event queues for these two types of 
events (not shown). In that embodiment, the program performs time sharing operations 
in which a predetermined percentage of the simulation time is used to process events 
from each queue. 

Other embodiments of simulation programs rank frequency channels on a per 
base station basis rather than on a per angular sector basis. 

Other embodiments of the invention will be apparent to those skilled in the art in 
light of the specification, drawings, and claims of this application. 
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